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A recently reported symmetry breaking of density profiles of fluid argon confined by two parallel 
solid walls of carbon dioxide is studied. The calculations are performed in the framework of a 
nonlocal density functional theory. It is shown that the existence of such asymmetrical solutions 
is restricted to a special choice for the adsorption potential, where the attraction of the solid-fluid 
interaction is reduced by the introduction of a hard-wall repulsion. The behavior as a function of 
the slit's width is also discussed. All the results are placed in the context of the current knowledge 
on this matter. 
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I. INTRODUCTION 

In a quite recent paper, Berim and Ruckenstein pj 
have reported symmetry breaking of the density profile 
of fluid argon (Ar) confined in a planar slit with identical 
walls of carbon dioxide (CO2). These authors claimed 
that a completely symmetric integral equation provides 
an asymmetric profile which has a lower free energy than 
that of the lowest symmetric solution leading to a sym- 
metry breaking phenomenon. It was assumed that the 
Ar atoms interact via a standard Lennard- Jones (LJ) po- 
tential characterized by the strength s ff and the atomic 
diameter aff. The presented results were obtained from 
calculations carried out with the smoothed density ap- 
proximation (SDA) version 0, [1] of the nonlocal density 
fmictional theory (DFT) in the case of a closed planar 
slit with an effective width of 15 <7//. The symmetry 
breaking was found for temperatures between the exper- 
imental triple point for Ar, Tt = 83.8 K, and a critical 
value Tsb = 106 K. At each temperature, it was deter- 
mined a range of average densities psbi < Pav < Psb2 
where the symmetry breaking occurs, outside this range 
a symmetric profile has the lowest free energy. 

As a matter of fact, the adsorption of fluid Ar on a 
solid substrate of CO2 was intensively studied for sev- 
eral decades (see, e.g., Refs. U andlSj). In 1977 Ebner 
and Saam analyzed phase transitions by assuming 
that atoms of the fluid interact with the solid wall via 
a 9-3 van der Walls potential (from here on denoted 
as ES potential) obtained from the assumption that Ar 
atoms interact with CO2 molecules via a LJ interac- 
tion with parameters Ssf and (Tsf- After this pioneer- 
ing work, a large amount of work has been devoted 
to study this system with different numerical and ana- 



lytical techniques. The attention was focused to ana- 
lyze features like: the oscillatory behavior of the den- 
sity profile which leads to a layered structure in the 
neighborhood of the flat substrate; the thin- to thick- 
film transitions; wetting properties and prewetting jumps 
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Berim and Ruckenstein have also adopted the ES po- 
tential, however, a hard-wall repulsion was introduced in 
their calculations. In practice, such a hard wall dimin- 
ishes the strength of the solid-fluid interaction. 

The investigation of symmetry breaking in physical 
systems is a very exciting issue. This is due to the fact 
that such a feature may have fundamental theoretical 
implications. In many fields of physics the discovery of 
symmetry breaking lead to significant advances in the 
theory. Therefore it is important to place the results of 
Ref. [1[ in the context of the current knowledge about 
adsorption into planar slits. 

Asymmetric solutions for fluids conflned in slits have 
been previously reported in the literature. About 20 
years ago a Dutch Collaboration has carried out cal- 
culations on the Delft Molecular Dynamics Processor 
(DMDP), which was specially designed for Molecular Dy- 
namics (MD) simulations of simple fluids [H, |2^. The 
results werepublished in a series of papers by Sikkenk 
et al. [H, Hi] and Nijmeijer et al. HI]. The sim- 
ulations were performed for a canonical ensemble with 
two types of particles, 2904 of one type for building a 
solid substrate and several thousand of the other type 
for composing the fluid adsorbate. The temperature of 
the system was kept a.t T* — ks T / e f f — 0.9 which is in 
between the fluid's triple-point temperature Tf ~ 0.7 and 
the critical temperature T* ~ 1.3. The width of the slit 
was taken as L = 29.1(7//, supposing that this distance 
be enough large to avoid any capillary effect. Such a sys- 
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tem can support solid-liquid (SL), solid-vapor(SV), and 
liquid- vapor (LV) interfaces. These authors have stud- 
ied wetting at LV coexistence by varying over a wide 
range the relative strength of the solid-fluid and fluid- 
fluid interactions defined by the ratio Er = £sf / £ff of 
the LJ parameters. The length scale of this interaction 
was taken as <Tsf = 0.941 ct//. For increasing from 
Er ~ 0.1 towards — 1-0 three cases were observed: 

(i) at low Er, symmetric profiles consisting of two SV 
interfaces and two LV interfaces are obtained, this sit- 
uation corresponds to a complete wall drying as can be 
seen in Fig. 1 of [111; 

(ii) at intermediate Er, asymmetric profiles consisting 
of a SL, a LV, and a SV interface are obtained, here the 
wall attraction is sufficiently strong to produce a partial 
wetting, i.e., to support a rather thick film on one wall 
while a SV interface is present near the other wall, this 
feature is shown in Fig. 2 of [28'|; 

(iii) for the largest Er, symmetric profiles consisting of 
two SL interfaces and two LV interfaces are obtained, 
now the strength is enough to wet both walls as can be 
seen in Fig. 3 of 28]. 

The structure of the profiles mentioned above depends 
on the balance of the involved surface tensions 7sl, Ilv, 
and which are related by the Young's law (see, e.g., 
Eq. (2.1) in Ref. 0) 

ISV = 1SL+1LV C0S9 , (1.1) 

where 9 is the contact angle. The latter quantity is de- 
fined as the angle between the wall and the interface be- 
tween the liquid and the vapor (see Fig. 1 in Ref. [29j). 
The transition from (i) to (ii) takes place at the drying 
point = 7:, whereas the transition from (ii) to (iii) takes 
place at the wetting point = 0. 

It is worth of notice that Velasco and Tarazona [s^l 
have carried out calculations in the framework of the SDA 
obtaining density profiles with the same structure to that 
reported in Refs. [25, 26, 27, 28]. The reader may look at 
Ref. [3l| for a further comparison between MD and DFT 
results. 

It is the aim of the present work to acquire a more 
accurate picture of the symmetry breaking reported in 
Ref. [l|. In so doing, we explore size effects by com- 
paring the results obtained for shts of widths 15 cr// 
and 30(7//. Next, we investigate the existence of sta- 
ble asymmetric solutions for the density profiles when 
the position of the hard-wall repulsion introduced in Ref. 
[l| is changed. When the location of this hard wall is 
moved, the strength of the adsorption potential is varied 
allowing a connection to the studies described in Refs. 
di, HEM m, im . Several properties of the obtained 
solutions are discussed. 

The paper is organized as follows. In Sec. |TT]we pro- 
vide a summary of the formalism underlying the present 
calculations. Special attention is devoted to the adsorp- 
tion potential. The results and its analysis are reported 
in Sec. IIIII Final remarks are given in Sec. II VI 



II. THE MODEL 

The properties of a fluid adsorbed by an inert solid 
substrate may be studied by analyzing the grand free 
energy |32] 

Vt = F-iiN, (2.1) 

where F is the Helmholtz free energy, /i the chemical 
potential, and N the number of particles of the adsorbate 

N = j p{r)dr . (2.2) 

Quantity F contains the energy due to the interaction 
between fluid atoms as well as the energy provided by the 
confining potential. In a DFT it is expressed in terms of 
the density profile p{r) 

F[p{r)] = F,^Mr)] + J dv p{v) Usf{r) . (2.3) 

Here J]nt [p(r)] is the intrinsic Helmholtz free energy func- 
tional and Usf{r) is the external potential produced by 
the slit's walls. 

This formulation is usually applied to systems de- 
scribed by the grand canonical ensemble, i.e., at con- 
stant volume V, temperature T, and chemical potential 
p. Such a situation corresponds to an open system in 
contact with reservoir which fixes T and p. A minimiza- 
tion of fl with respect to p{r) leads to the Euler-Lagrange 
equation for the density profile and the number of parti- 
cles may be evaluated with Eq. (|2.2p . For a closed system, 
i.e., a canonical ensemble with fixed N, one should treat 
p as an unknown Lagrange multiplier to be determined 
from the minimization procedure. 

A. Density functional theory 

Let us now summarize the DFT adopted for Fint[p{r)]. 
In the case of inhomogeneous classical fluids at tempera- 
ture T the intrinsic free energy functional is decomposed 
into two kind of contributions: 

(i) the ideal gas term i^id[jo(r)], which is given by the 
exact expression 

FiMr)] ^ ^bT j dv p{v) ha[{p)] 

= keT [ drp(r){ln[A3p(r)] - 1} , (2.4) 



with A = \J2tt It? /m ks T being the thermal de Broglie 
wavelength of a molecule of mass to; 

(ii) the excess term Fc-x[p{t)], which accounts for the 
interparticle interactions is a unique but unknown func- 
tional of the local density. For fluids with attractive in- 
teractions as the Lennard- Jones (LJ) one, the free energy 
is decomposed into the repulsive and attractive contribu- 
tions. The repulsive interactions are then approximated 
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by a hard-sphere functional with a certain choice of the 
hard-sphere diameter dfjs 



^Hs[p(r)] = / drp(r)/Hs[/3(r);dHs] 



(2.5) 



whereas the attractive interactions are treated in most 
cases in a mean field fashion 

i^attr[p(r)] ^\ J J dvdvlp{v) p{vl) $attr(| ^ - Yl \) . 

(2.6) 

Here $attr(^ =| v — ri |) is the attractive part of the LJ 
potential. 

In summary, the intrinsic Helmholtz free energy func- 
tional may be expressed as 

Fi„t[p(r)] - i^id[p(r)]+Fe,[p(r)] 

= i^id[p(r)] + FHs[p(r)] + Fatt.[p(r)] (2.7) 

The free energy functional for hard spheres plays a cen- 
tral role in DFT. Expressions for /Hs[p(r); ^hs] may be 
taken from the Percus-Yevick [s^ or Carnahan-Starling 
(CS) |3J] approximations for the equation of state (EOS) 
of a uniform non-attractive hard-sphere fluid (see, e.g. 
Ref. [1^). In a nonlocal DFT this quantity is evaluated 
as a function of a conveniently averaged density p{v). 

For the calculations performed in the present work we 
used the same SDA formalism adopted in the paper of 
Berim and Ruckenstein [l[. In this approach developed 
by Tarazona 0, [1] , the excess of free energy density of 
hard spheres is written according to the semi-empirical 
quasi-exact CS expression 



/Hs[p(r); dns]] = ./csW = ksT 



??(4-3^) 



(2.8) 



Here fj — p{r) Vhs is the packing fraction, where the fac- 
tor Vhs — ■"'c^Hs/^ volume of a hard sphere. The 
smoothed density p(r) is defined as 



p(r) = J dr p(r) w[| r — r/ |;p(r) 



(2.9) 



with the following weighting function: 



w[\ r — r/ I; /c(r)] 



wo[\ r~r/\]+wi[\r-r/ \] p(r) 
+W2[\r-r/\]p'^{r) . (2.10) 



The expansion coefficients woir), wi{r), and W2{r) are 
density independent and its expressions as a function of 
r = I r — r/ I are given in the Appendix of Ref. i3|. 

To account for the fluid-fluid interaction we adopted, 
as in Ref. 1], the spherically symmetric L-J potential 
given in Eq. (2) of Ref. 8 



4'attr(?') 



4e// 




^^"^'^//(2..11) 
if r < aff 



where cr^/ is the hard-core diameter of the fluid. The au- 
thors of Ref. 8 have used this L-J version just for studying 



the adsorption of Ar on CO2 and it has been also utilized 
in several subsequent works on this system. The values of 
the interaction parameters for Ar are Eff/ks = 119.76 K 
and <Tff = 3.405 A. 



B. The Euler-Lagrange equation 

The equilibrium density profile p(r) of the fluid ad- 
sorbed in a closed slit is determined by a minimization 
of the free energy with respect to density variations with 
the constraint of a fixed number of particles N 



5p{r) 



Fi„t[p(r)]+ / drip{vt)[Usf{vi)-p] 







(2.12) 



Here, i.e., for an ensemble with fixed V , T, and N, the La- 
grange multiplier fi is an unknown quantity which should 
be determined from the constraint. It plays a role of a 
chemical potential but off the liquid-vapor coexistence 
conditions. Hence, it is not necessarily equal to fJ-coex of 
an open slit in equilibrium with a reservoir at tempera- 
ture T (see, e.g., Ref. IH). 

In the case of a planar symmetry where the flat walls 
exhibit an inflnite extent in the x and y directions the 
proflle depends only of the coordinate z perpendicular 
to the substrate. For this geometry the variation of Eq. 
(|2.12p yields the following Euler-Lagrange (EL) equation 



Spiz) 



Sp{z) 

+ J dzlp{zl) iattrd Z- Zl\) + U,f{z)^ p . 

(2.13) 

For a slit of effective width tyj this EL equation may be 
cast into the form 



kBTlnlk" p{z)] + Q{z)= p , 



(2.14) 



where 



Q{z) = ksT 



Afj{z) - ■iff{z) 



^i. rj^^dls , , 4-2r7(z/) SpizQ 

+ / dzfp{z/) i-attrd Z~Z/\) + U,f{z) (2.15) 

Jo 

The number of particles per unit area of one wall of the 
slit is 



p{z) dz . 



(2.16) 



In order to get solutions for p{z) it is useful to rewrite 
Eq. ((2T4)) as 



p{z) = po exp - 



Q{z) 
ksT 



(2.17) 
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with 



Po 



1 

A3 



exp 



(2.18) 



The relation between ^ and Ns is obtained by substitut- 
ing Eq. (|2.17p into the constraint given by Eq. (|2.16p 



fcs T In 



1 



dz exp — 



(2.19) 



For the calculations carried out in the present work we 
set dns = <^ff as it was done in Ref. 111. 

The asymmetry of the density profiles is measured by 
the parameter 



2iV., 



dz I p{z) - p{luj - z) I 



(2.20) 



According to this definition, if the profile is completely 
asymmetrical about the middle of the slit [p{z < Iw/'^) 
and p{z > 1^/2) = 0] this parameter becomes unity, 
while for symmetric solutions it vanishes. 
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FIG. 1: Adsorption potential as a function of the distance 
from the real wall. The solid curve is the potential of Berim 
and Ruckenstein close to the left wall given by Eq. (|2.23|l . 
while the dashed curve is that of Nilson and Griffiths given 
by Eq. (p:25)l . 



C. Adsorption potential 

The model van der Waals (9-3) potential proposed by 
Ebner and Saam i.e. the ES potential, is 



Usfiz) 



2lT ^ 





if z > 0, 



if z < 



'(2.21) 



with Eef f — Esf Ps cr^f being the effective strength, was 
adopted for almost all the abovementioned studies of the 
adsorption of Ar atoms on a fiat wall of solid CO2. The 
exception is the experimental and theoretical investiga- 
tion performed by Mistura et al. (2^ . where a more re- 
alistic adsorption potential calculated on the basis an 
ah initio expansion of Marshall et al. [st'I was used. 
The ES expression is obtained when one assumes that 
Ar atoms interact with CO2 atoms via a Lennard-Jones 
(12-6) potential and subsequently integrates this poten- 
tial over a continuum of CO2 substrate atoms with a re- 
duced density p* = Ps<^sf = 0.988. The cross-parameters 
of the potential are determined by using the Lorentz- 
Bcrthclot rules. So that, the van der Waals strength Esf 
is the square root of the product of the argon and CO2 
van der Waals strengths, while the hard-core diameter 
(Ts/ is the mean of the argon and CO2 hard-core diam- 
eters, while The parameters evaluated in this way are 
Esf/kB = 153 K and asf = 3.727 A. 

Berim and Ruckenstein [l| have investigated the Ar- 
CO2 system utilizing, in principle, the ES potential. 
However, by looking at their paper one realizes that ac- 
cording to Eq. (A5) of the Appendix 



Usfl{z) = Y^e// 



2 

15 



z + asf 



(7sf 



z + asf 



(2.22) 



which accounts for the solid-fiuid interaction at one of 
the walls, a hard-wall repulsion was located at a distance 
asf from the real wall of the slit. In agreement with this 
assumption, the total confining potential exerted on Ar 
atoms by the two walls separated by a distance L was 
expressed as 

Usf{z) = Usfi{z + asf) + Usf2{(-w -z + asf) ■ (2.23) 

Here the effective width of the slit is 

i^^L^2asf. (2.24) 

This scenario is depicted in Fig. 2 of Ref. T. In this con- 
text, it is interesting to notice that Nilson and Griffiths 
[3^ in order to study the adsorption of a fiuid in a pla- 
nar slit have written in their Eq. (10) the total fiuid-solid 
potential as 



Usf{z) 



Usfi{z+^) 



Usf2{iw ~ z + ^) , (2.25) 



i.e., locating a hard- wall repulsion at a distance asf/2 
from the substrate. In this case, the effective width is 



(Jsf 



(2.26) 



as it is shown in Fig. 2 of Ref. l38l . 

In Fig. [1] we compare the potentials outlined in the 
previous paragraph. The comparison is restricted to the 
region close to the substrate. The quantity ^ is the per- 
pendicular distance from the real wall being 



z -f- asf for Berim — Ruckenstein 



z + 



for Nilson — Griffiths. 



(2.27) 



One may realize that Eq. (|2.25p retains the "soft" re- 
pulsion \Usf{z) oc ((Tsz/z)^], while Eq. (|2.23p cuts the 
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TABLE L Values of the Helmholtz free energies Fsym, Fasym, 
and Fcap of the symmetric, asymmetric and capillary solu- 
tions, respectively, obtained with v = 1 for the slit £^ = 15 
at T = 87 K. The free energies are given for several average 

2 

//■ 



densities in the range p*t,i < Pav < plb2 in units of ks T/a 





Fsym 




Fasym 


Fcap 


Pav 


BR" 


PW'' 


BR° 


PW'' 


PW*" 


0.1546 


-26.59 


-26.62 


-26.67 


-26.70 




0.2319 


-39.66 


-39.69 


-39.86 


-39.88 


-39.12 


0.3092 


-52.81 


-52.85 


-53.06 


-53.08 


-52.32 


0.3865 


-66.01 


-66.04 


-66.26 


-66.29 


-65.65 


0.4638 


-79.22 


-79.28 


-79.48 


-79.51 


-79.26 



" Data taken from jj]. 
Calculated in the present work. 



III. NUMERICAL RESULTS AND DISCUSSION 

Let us now describe the obtained results. The EL equa- 
tion (I2.14P was solved at fixed and T for a given num- 
ber of particles per unit area Ng. The latter quantity de- 
termines an average fluid density pav = Ng/iw A widely 
used computational algorithm consisting of a numerical 
iteration of the coupled Eqs. (|2.17p - (|2.19p was applied. 
This procedure yields the density profile p{z) and the 
value of the Lagrange multiplier /i. The convergence of 
the solutions are measured by the difference between two 
consecutive profiles 
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'ff 



dz 



1 2 



Pi+liz) - p,{z) 



(3.1) 
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where pi(z^ is the density profile after the i-th iteration, 
and by the quantity 
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FIG. 2: Lagrange multiplier /i as a function of average density. 
The solid curve are results for symmetric film solutions. The 
dashed curve stands for values of asymmetric film solutions 
which occur in the range p^sh\ — Pav ^ P*sb2 ■ The dot-dashed 
curve corresponds to drying-CC like solutions. 



potential before the minimum be reached. This feature 
produces important effects on the behavior of the density 
profiles. In fact, the calculations performed by Berim and 
Ruckenstein [l| yielded density profiles with p{z — 0) and 
p{z — £m) different from zero indicating that the fluid is 
in contact with the hard walls, while in the case of Nilson 
and Griffiths [s^ the fluid forms a well defined first layer 
separated from the wall. 

In the present work we shall analyze the evolution of 
asymmetric solutions when the total adsorption potential 
is written as 



Ugfiz) - Usfiiz + ^) + C/./2(C - z + ^) , (2.28) 



and the parameter v varies from 1 to 2. In doing so, one 
goes from Eq. (|2.23p towards Eq. (|2.25p increasing the 
strength of the solid-fiuid attraction. 



(52 = 1 - — 



dz p(z) , 



(3.2) 



accounting for the deviation from the required Ns- 

In practice, for the calculations it is convenient to use 
dimensionless variables: z* = z/uff for the distance, 
p* — pcr^ff for the densities, and T* = kgT/eff for the 
temperature. In these units the average density becomes 
p*^ = Nsajf/il = N*/£l. For the numerical task, the 
region of integration [0, was divided into a grid of 
equal intervals Az* = 0.02, i.e., a grid with 50 points per 
atomic diameter Uff. It is worthwhile to notice that in 
the work of Berim and Ruckenstein the number of grid 
points was taken equal to 10 per atomic diameters. If the 
obtained profile did not change with increasing precision 
from Si « 10~® to Si « 10~^^, then it was accepted as a 
solution of the coupled integral equations. 

In a first step, we studied the same systems treated in 
detail by Berim and Ruckenstein Hence, we set ly = 1 
and solved the EL equation for a slit with an effective 
width C = 15 at T = 87 K (T* = 0.73) for a series 
of average fluid density p*^ — iV*/^^. The ground state 
solutions yield symmetric density profiles for p*„ < 0.1 = 
P*sbi and for p*„ > 0.514 = p*^2, while in the range < 
Pav — P*sb2 ^^'^ ground-state solutions provide asymmetric 
density profiles. This is due to the fact that in such a 
regime the asymmetric solutions have lower free energy 
than the symmetric ones. The free energies calculated for 
some selected /O*^ are listed in Table [J together with the 
results obtained by Berim and Ruckenstein py] . A glance 
at this table indicates a good agreement between both 
sets of values. In order to get symmetric solutions in the 
range 0.1 < p*„ < 0.513 one must explicitly impose such 
a condition to Eqs. ((2l7|) - ([2T9l) . 

The Lagrange multiplier p (equivalent to the chemical 
potential in the case of open slits) is displayed in Fig. [2] 
as a function of average density. We show the results for 
a wider range of p*^ than it is done in Fig. 9 of Ref. [H. 
Figure [2] clearly indicates that the asymmetric solutions 
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2.5 5 7.5 10 12.5 15 



FIG. 3: Solutions of the EL equation (|2.14p for the average 
density p*„ — 0.3865. The solid curve corresponds to the 
lowest symmetric solution, the dashed to the asymmetric one, 
and the dash-dotted to the drying-CC like one. 

1.25| ■ ■ • ■ ■ 1 
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z/a„ 

FIG. 4: Asymmetric density profiles as a function of the dis- 
tance from the hard-wall repulsion. It is shown the evo- 
lution for increasing density average at u — 1. Data for 
= 0.1063, 0.1159, 0.1546, 0.1932, 0.2319, 0.3092, 0.3892, 
0.4638, 0.5135, and 0.5139 are displayed. These values cover 
the whole range p*i,i < pav < P*sb2 where there are asymmetric 
solutions. 

occur in the domain where /i is degenerate, namely, the 
same value of ^ corresponds to different p*„. That is just 
the regime where in the case of an open slit the equal- 
area Maxwell construction should be applied in order to 
determine the chemical potential (cf., e.g., Fig. 2 in Ref. 
HB)- In addition, it is worthwhile to notice that at p*fj2 
there is an abrupt jump of /i, we shall come back to this 
feature below. 

If one manages conveniently the EL equation it is pos- 
sible to get another kind of symmetric solutions in some 
region of The free energy of such new solutions is 
included in Table U and the corresponding multiplier fi 
is plotted in Fig. [21 Figure [3] shows the spacial distribu- 
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FIG. 5: Asymmetry parameter defined in Eq. (|2.20p as a func- 
tion of average density. The successive curves are results ob- 
tained for slits with different effective widths labeled by 
For each width the asymmetric solutions occur in different 
ranges p*6i < pL < p*62- 

tion of all three states listed in Table [J for p* ,j = 0.3865. 
A direct comparison indicates that these three profiles 
correspond to the drying^ the one-wall-wetting , and the 
two-walls-wetting cases displayed, respectively, in Figs. 
1-3 of Ref. 2^1 . From a glance at Table [T] one realizes 
that as long as the film wetting solutions exist the drying 
one is a symmetric excited state. In the regime p*^ > p*^2 
the drying profile shown in Fig. [3] becomes the capillary 
condensation (CC) solution. 

Figure 2] shows a series of asymmetric density profiles 
p*{z) = p{z)a^j:. This sequence indicates that for in- 
creasing average density, starting from the profile de- 
noted as 1 (p*„ = 0.1063), the number of oscillations near 
the left wall as well as its amplitudes increase. This trend 
continues up to the profile 2 (p*^ = 0.4638), where the 
peaks of the oscillating structure begin to decrease. Fur- 
thermore, profile 3 corresponds to the biggest asymmetric 
solution, for larger p*„ there is a jump to the symmetric 
profile 4. The latter behavior indicates a transition to the 
so-called CC phase, i.e., a transition to a situation where 
the slit is full of liquid argon. The jump of ji addressed 
in the previous paragraph is also a manifestation of this 
thick-film to CC transition. 

The asymmetry of the profiles for v — I displayed in 
Fig. |4] was measured by the parameter A a? introduced 
in Eq. (|2.20p . The results are shown in Fig. [5] These 
values are essentially the same as that of the equivalent 
curve displayed in Fig. 5 of Ref. |l|. As mentioned in 
the introduction, in Refs. [H, [H, [23, [H, H [U it is 
emphasized that a slit of width = 30 is appropriate 
to study wetting because it is large enough to avoid any 
confinement effect. Therefore, it becomes of interest to 
solve the EL equation for such a big slit and to compare 
the results with that obtained for = 15. The asymme- 
try parameter evaluated for the £^ = 30 slit at T = 87 
K is also plotted in Fig. O A comparison of the curves 
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FIG. 6: Asymmetry parameter of density profiles, defined 
in Eq. (|2.20[) . as a function of average density. Going from 
tfie outer curve towards the inner one the data correspond, 
respectively, to u ^ 1, 1.136, 1.156 and 1.170. Uu> 1.18 for 
all Pod one gets Ajv = . 



I p =0.1932 



0.6 




T = 87 K 
I =15o„ I 

W ff * 



J t 



7.5 10 12.5 15 



FIG. 8: Evolution of the asymmetric solution for the average 
density p*„ = 0.1932 in terms of the factor u. The solid curve 
corresponds to = 1, the dashed to 1.136, the dash-dotted 
to 1.156, and the dotted to 1.2. 





T = 87K 




1 =15 0., 

w ff 


p* =0.1932 





1 1.04 1.08 1.12 1.16 1.2 



V 

FIG. 7: Asymmetry parameter of density profiles, defined in 
Eq. p.20|l . obtained for the average density p^D = 0.1932 as 
a function of the factor v. 



labeled by 30 and 15 indicates that its shapes do not 
differ significantly, in particular, in the regime of large 
coverage 0.35 < p*^ < 0.51 the values of Ajv are almost 
equal. Hence, we can state that the capillary effects in 
the case of the moderately thick slit of = 15 are not 
very important. 

For the sake of completeness we also analyzed slits of 
£^ < 15atr = 87K. The asymmetry parameters evalu- 
ated for €^ = 12, 7.5 and 6 are included in Fig. (5] It is 
clear that these curves exhibit everywhere an increasing 
departure from the — 30 results. Finally, for = 5.5 
the parameter A^v becomes zero for all p*^ indicating 
that the asymmetric solution disappears. 

In the next step, we analyzed the change of the prop- 
erties described above when the parameter v is taken 
larger than unity. The variation of v is performed in 



such a way that the effective width of the slit is always 
kept at = 15, hence, only the potential is slightly 
changed mainly near the walls. At this point the tem- 
perature was still kept at T = 87 K. It was found that 
for increasing values of v the range of average density 
Psbi ^ P*av — P*sb2 where there are symmetry breaking 
decreases. This effect might be observed in successive 
plots of /xvs.p*„, however, we prefer to report directly 
the evolution of the parameter Ajv . Figure [6] shows how 
this parameter decreases with increasing ly. From this fig- 
ure one may conclude that symmetry breaking persists 
at most for p*^ « 0.2. 

Figure [7] shows the asymmetry parameter A^r as a 
function of the factor for the average density p*,y — 
0.1932. These data indicate that the asymmetric solu- 
tion already disappears for a critical value — 1.18. 
The evolution of the density profiles from asymmetric to 
symmetric species is displayed if Fig. [51 In this draw- 
ing one may observe the diffusion of fluid argon from 
the neighborhood of the left wall towards the right one. 
This process continues until a symmetric density profile 
is formed for ~ 1.2. The described evolution of p*{z) 
is determined by the behavior of the surface tensions. If 
one uses the Young's law given in Eq. Ijl.ip the total sur- 
face excess energy of asymmetric profiles may be written 
as 

Itot = JSL+lLV+lSV = 2 75L + 7Ly (l + cos^) , (3.3) 

with cos 6* = {"^sv ~1sl) IlLV < 1- By increasing enough 
the attraction the equality ^ysv — ISL = ILV is reached 
yielding cos^ — 1, then the system undergoes to a tran- 
sition to a symmetric profile with 

7tot = 2 75i + 2 7Ly . (3.4) 

In this case both walls of the slit are wet. 
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Let us now look if something special is going on for 
the adsorption potential at ~ 1.18. As a matter of 
fact, it becomes important to explore where the hard- 
wall repulsion introduced in Rcf. 1 is located for such a 
critical value. The position of the left hard-wall repulsion 
with respect to the real wall according to Fig. [1] is given 
by 



HW 



(3.5) 



For a slit of effective width £^ — 15 the minimum of the 
total adsorption potential given by Eq. (|2.25p is to a very 
good approximation determined by the minimum of the 
ES potential exerted by the left wall. It is located at 



1/6 



The ratio of these quantities becomes 



HW 



1/6 



1.013 



(3.6) 



(3.7) 



Since this ratio is larger than unity, it indicates that for 
Vc the minimum of the ES potential is reached inside the 
slit. 

We also analyzed the occurrence of asymmetric solu- 
tions for 1 < z/ < 2 at temperatures T > 87 K. It was also 
found that the symmetry breaking disappears at some 
critical value of v. Such a behavior may be expected if 
one takes into account that for temperatures larger than 
T = 87 K even for v = \ the range psbi < Pav < Psb2 is 
smaller (see Fig. 5 of Ref. [l]) . 

It is worthwhile to notice that a reliable ab-initio po- 
tential utilized by Mistura et al. [l^] for investigating 
the adsorption of Ar on CO2 exhibits an even stronger 
attraction than the ES p otential. This feature can be 
observed in Fig. 3 of [24I. Hence, one would not expect 
any symmetry breaking in the case of such a realistic 
adsorption potential. 



IV. FINAL REMARKS 

We reexamined the symmetry breaking found very re- 
cently by Berim and Ruckenstein [3| in a study of the 
adsorption of argon in a closed slit with identical walls of 
carbon dioxide. It is important to stress that these au- 
thors have introduced hard-wall repulsions at distances 
dsR = o'sf from the real walls of the slit as shown in Fig. 
[T] (see also Fig. 2 in Ref. [IJ), reducing in such a way 
the strength of the adopted ES potential. Stable asym- 
metric solutions were obtained in the case of a slit with 
effective width £^ = 15 for temperatures in the range 
87 < T < 106 K. 

The present study was mainly devoted to establish how 
robust are the asymmetric solutions against changes of 
the adsorption potential. However, in addition, the be- 
havior of the asymmetry parameter was also evaluated 



for slits of different widths. The calculations have been 
carried out using the same nonlocal formalism as that 
adopted in Ref. [l|, namely, the SDA density functional 
theory proposed by Tarazona 

Since in a pioneering series of works a Dutch Collabo- 
ration have previously found the symmetry breaking by 
studying adsorption in a slit of width i^, — 30, see Refs. 
[25I . [27I . i28] , we performed a comparison of the asym- 
metry parameter Ajv evaluated for slits of £^ = 15 and 
30 at T = 87 K. The difference between the results is 
rather small. Further calculations showed that the effect 
due to confinement begins to be important for slits with 
C < 12. 

Focusing the analysis on the slit of = 15, in a rather 
complete plot of pvs. p*^ we show clearly the regime 
where the symmetry breaking occurs and, in addition, 
we also display the Lagrange multiplier p for a CC like 
solution. Furthermore, it is shown that the spacial shape 
of the obtained solutions correspond to the three sorts of 
profiles discussed in Refs. [H, iM [13, . The strength 
of the attraction determines which one of that profiles is 
the stable solution. 

By shifting the hard-wall repulsions introduced by 
Berim and Ruckenstein [l| towards the real walls the 
strength of the attraction is increased producing changes 
in the balance of surface tensions. It was found that at 
T = 87 K already for the critical distance c?c — ^Bfl/l-Z 
the asymmetric solutions disappear. At this tempera- 
ture, close to the triple point, the asymmetry parameter 
Ajv reached its largest values in the study of Ref. U. For 
higher temperatures the symmetry breaking disappears 
even more rapidly. 

Finally we can state that from the present study it 
is possible to conclude that the symmetry breaking re- 
ported in Ref. [l] can be understood in terms of findings 
described in Refs. [1^, [26, 27, 28]. By locating hard-wall 
repulsions the authors of Ref. 1 diminish the attraction of 
the ES adsorption potential exerted on the fluid causing 
the entrance of the system in a one- wall- wetting regime 
(asymmetric profiles) corresponding to Fig. 2 of Ref. [28j . 
By moving the hard- wall repulsions towards the real walls 
the attraction increases monotonically causing eventu- 
ally the entrance of the system in the two-walls-wetting 
regime with symmetric solution of the type displayed in 
Fig. 3 of Ref. HI. 

Furthermore, a reliable ab-initio potential for the in- 
teraction of fluid Ar with a structureless smooth wall of 
CO2 like that adopted by Mistura et al. is even 

stronger than the ES one. Therefore, it is possible to in- 
fer that no symmetry breaking would be expeted for real 
Ar/C02 systems. 
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